Optical detection of fractional particle number in an atomic Fermi-Dirac gas 
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We study theoretically a Fermi-Dirac atomic gas in a one-dimensional optical lattice coupled to a 
coherent electromagnetic field with a topologically nontrivial soliton phase profile. We argue that the 
resulting fractional eigenvalues of the particle number operator can be detected via light scattering. 
This could be a truly quantum mechanical measurement of the particle number fractionalization in 
a dilute atomic gas. 
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It has been known for a while that, in the presence 
of a topologically nontrivial bosonic background field, 
fermionic particles may carry a fractional part of an ele- 
mentary quantum number |3I ■ In the condensed mat- 
ter regime this phenomenon was introduced Q to de- 
scribe conjugated polymers. The existence of fractionally 
charged excitations in the polymers is typically demon- 
strated indirectly by detecting the reversed spin-charge 
relation The fractional quantum HaU effect (FQHE) 
can also be explained by invoking quasiparticles, each 
with a fraction of an electron's charge The fluctua- 
tions of the tunneling current in low-temperature FQHE 
regime have been measured Interpreting the cur- 
rent shot noise according to the Johnson-Nyquist formula 
duly suggests that the current is carried by the fractional 
Laughlin quasiparticles. Analogous experiments have de- 
termined the fractional expectation value of the charge 
in FQHE in the Coulomb blockade regime 

We have earlier proposed a system of Fermi-Dirac (FD) 
atoms in an optical lattice that should display fractional 
atom numbers In the present Letter we argue that 
usual optical methods such as phase contrast imaging and 
measurements of the intensity of light scattered by the 
atoms extract information about the fractional fermion 
number. The technical challenges are severe, but in prin- 
ciple both the expectation value and the fluctuations of 
the atom number are accessible to experiments. 

Briefly, we consider a FD atom |9| with a A scheme 
for two active states (say, Zeeman states in different hy- 
perfine levels) that can be coupled by one or two-photon 
electromagnetic (ern) transitions . The atoms reside in 
a ID optical lattice [l^ that holds the two states at alter- 
nating sites A/4 apart, where A denotes the wavelength 
of lattice light. By making use of the em transitions, we 
assume, it is possible to make the atoms hop between the 
adjacent sites so that they at the same time change their 
internal state. The lattice Hamiltonian is 



H 



E 



(1) 



and Kk are the are hopping matrix elements tailored to 
suit our purposes. The literature is replete with stud- 
ies of similar systems 0, , but our treatment is 
unusual in that we never resort to the continuum limit. 

It is easy to see that the matrix Ukp whose columns are 
the orthonormal eigenvectors of the eigenvalue problem 



(2) 



may be employed to diagonalize the Hamiltonian . In 
terms of the new fermion operators 7p = UkpCk, we 
have H/h = J^p ^piIip ■ Without a loss of generality, 
the couplings Kk are assumed real, and so we take U real 
and orthogonal; U^^ = Upk- In this paper all calculations 
are done directly numerically. 

Even though the optical lattice may be part of a larger 
trap and the wider trap could generate interesting physics 
in its own right, we simplify by putting 4 = 0. The 
fractional charge arises from certain types of defects in 
the couplings Kk- We illustrate by assuming a dimerized 
lattice generated by the coupling matrix element that al- 
ternates from site to site between two values a -f /i and 
a — fj,, except that at the center of the lattice there is a 
defect such that the same coupling matrix element ap- 
pears twice. We take the number of lattice sites to be 
Ns = 2Nh + 1 = An + 1, where n is an integer itself, and 
number the sites with integers ranging from —Nh to Nh- 
For illustration, pick n = 2, use an x to denote a lattice 
site, and ± the couplings a ± /i, then our lattice with the 
couplings reads 



X — X — X 



(3) 



Here Ck is the annihilation operator for the fermionic 
atoms at site fc, hSk is the energy of the atoms at fc, 



It is then easy to see from the structure of Eq. (0) that 
if Lu is an eigenvalue, then so is —lo; and the eigenvectors 
transform into one another by inverting the sign of ev- 
ery second component. We will label the eigenvectors as 
—Nh . . . Nh in ascending order of frequency, and assign 
the labels ±p to such ± pair of states. Correspondingly, 
the transformation matrix U satisfies \Ukp\ = \Uk,-p\- 

But under our assumptions, the number of eigenvalues 
and eigenstates is odd. The ± symmetry implies that an 
odd number of the eigenvalues must equal zero. Except 
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for special values of the couplings a and /z, there is one 
zero eigenvalue. We call the corresponding eigenstate 
the zero state. Provided a and have the same sign and 
\a\ > l/^l, all odd components in the zero state equal zero 
and the even components are of the form 



The zero state becomes the narrower, the closer in abso- 
lute value a and /i are. 

The dimerized optical lattice resulting from the alter- 
nating pattern of the hopping matrix elements causes the 
single-particle density of states to acquire an energy gap, 
which in the limit Ns ^ oo equals 4?i|/i| . The zero state 
is located at the center of the gap. The resulting exci- 
tations at half the gap energy could be detected by res- 
onance spectroscopy. This provides indirect evidence of 
fractionalization, as in the polymer systems 0] . Because 
in our scheme [3] the gap is proportional to the ampli- 
tude of the em field inducing the hopping, the size of the 
energy gap can be controlled experimentally. The effec- 
tive zero temperature limit, ^ ksT/h, might then be 
reached under a variety of experimental conditions. 

Suppose next that the system is at zero temperature, 
and contains Nf = Nh + 1 fcrmions. The exact eigen- 
states p are then filled up to zero state and empty at 
higher energies, with occupation numbers = or 1. 
The number operator for the fermions at site k corre- 
spondingly reads c\ck = Y.pq UkpUkqlllq, so the expec- 
tation value of the fermion number at the site k is 

Nh 
p=-Nh P=-Nh 

= \ + \\Uto?- (5) 

The second equality is based on the symmetry \Ukp\ ~ 
|[/fc._p|, and the third on the orthogonality of the ma- 
trix U . By virtue of the same orthogonality, localized 
with the zero state there is a lump with i \Uko\^ — \ 
fermions on top of a uniform background of half a fermion 
per site. This lump is the celebrated half of a fermion. 
So far we only deal with the expectation values of the 
atom numbers, but we will demonstrate shortly that the 
fluctuations in atom number can be small as well. 

Fractionilization is a more robust phenomenon than 
our discussion may let on. Something akin to a localized 
zero state occurs as soon as the regular alternation of 
the couplings between adjacent states gets out of rhythm 
around a defect. In particular, the defect does not have to 
be confined to one lattice site, which might make the ex- 
periments easier. The half-fermion is localized, so it does 
not critically depend on the number or parity of sites, 
and not even on the exact number of the fermions. We 
will enumerate such variations of the theme elsewhere. 

We next turn to the optical detection of the FD gas. 
We assume that far off-resonant light excites the atoms, 



whereupon the ID optical lattice may be considered op- 
tically thin. We take the sources residing at each site to 
be much smaller than the wavelength of light. Then the 
(positive frequency part of the) field operator for scat- 
tered light is Q E'^ = CJ2k '^kc\ck, where C is a con- 
stant containing the overall intensity scale of the driving 
light. Henceforth we scale so that C — \. The factors ak 
include aspects such as intensity and phase profile of the 
driving light, effects of the spin state at each site k on 
the light-atom coupling, and propagation phases of light 
from the lattice site to the point of observation. 

In forward scattering and variations thereof such as 
phase contrast imaging, the scattered and the incoming 
light interfere. The ultimate measurement of the inten- 
sity in effect records the expectation value of the electric 
field E = (E^)- The observable at the detector is 

E = J2'^k{clck) =J2^i'Ukpr^p- (6) 

k kp 

This is a linear combination of the expectation values of 
the numbers of fermions at each lattice site with the co- 
efficients ak , which are to some extent under the control 
of the experimenter. In the absence of interference with 
the incoming light, the simplest observable is light inten- 
sity / = {E^E^). The detector then probes the quantity 
characteristic of the FD statistic 

/ = E"^"'(4c4c;)-|i?l' + A/; (7) 

kl 

AI = ^ alaiUkpUkqUipUiqnp{l - Uq) . (8) 

klpq 

We now construct a numerical example about the use 
of forward scattering to detect the fractional particle. We 
make use of the fact that the fermion species at the alter- 
nating lattice sites are likely to be different. We assume 
that a given driving light is far above resonance with the 
fluorescing transition in one species and far below res- 
onance in the other, and that the matrix elements for 
dipole transitions are comparable. It is then possible to 
find a tuning of the laser such that the intensity of scat- 
tered light is the same for both species. Moreover, the 
light scattered by the two species are out of phase by tt, 
and out of phase with the incident light by tt/2. With the 
usual tricks of phase contrast imaging, the relative phase 
of incident and scattered light is then adjusted so that 
that in interference light from one species directly adds 
to the incident light, and light from the other species 
subtracts. 

The second element of the argument is a rudimentary 
model for an imaging system with a finite aperture. Let 
us assume the geometry has been arranged in such a 
way that all Fourier components of light in the plane 
of the aperture up to the absolute value K are passed, 
the rest blocked. The effect on imaging from the ob- 
ject plane to image plane can be analyzed by taking the 
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FIG. 1: Optical image (dashed line) of a zero state carrying 
half of a fermion (solid line) for a specific imaging system as 
discussed in the text. The size of the soliton is set by the 
choices Ns = 129, Nf = 65, = 0.1 a. 

Fourier transform of the object, filtering with the multi- 
phcr d{K — |kj_|), and transforming back. The filtered 
image of the lattice filtered is then proportional to 

We choose the parameters n = 0.1a, and the numbers 
of sites and fermions Ng = 129 and Nj — 65. We take the 
numerical aperture F ^ 1 for the imaging system, and 
the corresponding maximum possible cutoff wave number 
K = 27r/(\/5A). In Fig. [H we plot the optically imaged 
fermion lattice along the line of the atoms (dashed line), 
and the number of fcrmionic atoms in excess of the av- 
erage occupation number i for the even-numbered sites 
that carry the zero state (solid line), as obtained from 
Eq. |SJ). The curves are normalized so that the max- 
ima overlap. The imaging system picks up a resolution 
rounded version the half- fermion hump. 

In fact, phase contrast imaging has been used for 
nondestructive monitoring of a Bose-Einstein conden- 
sate 01 , and the absorption of a single trapped ion has 
been detected experimentally While a lot of as- 

sumptions went into our specific example, an experiment 
along these lines should be feasible with the technology 
available today. 

With illumination of the optical lattice by a focused 
light beam and detection of scattered intensity in a di- 
rection of constructive interference, it is in principle also 
possible to realize a situation in which the weights ap- 
proximately make a Gaussian distribution around the 
zero state, = e"^'^/™-' . In such a case the observable 
TV = is just a linear combination of the occupation 
numbers of the lattice sites, the quantity E is the expec- 
tation value thereof, and AJ is nothing but the squared 
fluctuations of N, AJ = (AN)^ = (N^) - (N)^. 

To illustrate, we take a lattice with Ng = 1025 sites, 
pick the parameters fi = 0.1 a, put in iV^ = 513 fermions 
so that the zero state is the last filled stated, and find 
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FIG. 2: Fermion number fluctuations AA^ under a Gaussian 
envelope of width w, for Ns = 1025, Nf = 65, /i = 0.1 a. 
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FIG. 3: The intensity of light scattered from the optical lattice 
if the fermion numbers did not fluctuate {\E\'^), and the addi- 
tional intensity due to fermion number fluctuations ((AA'^)'^), 
as a function of the size of the focus w of the driving light 
given in the lattice units of A/4. We also display the added 
intensity {AN)^ that would result if fermion number fluctu- 
ations were uncorrelated between adjacent lattice sites. The 
soliton parameters are Ns = 129, Nf = 65, fi = 0.9 a. 

the rms fluctuation of the fermion number AA'^ as a func- 
tion of the width of the weight function w. The result 
is shown in Fig. |21on a log-log plot. The notch around 
w = 1 indicates that at this point the weight factors 
start to cover several lattice sites. Another break in the 
curve is seen at about w = 10, when the weight function 
covers the whole zero state. Thereafter the fluctuations 
behave as A7V oc The fermion number N un- 

der the weight function becomes more sharply defined as 
the region for averaging grows broader. Finally, at w of 
a few hundred, the weights effectively cover the en- 
tire lattice. The fluctuations then decrease even faster 
with increasing w, as is appropriate for the fixed fermion 
number in the lattice as a whole. 

In the standard half-integer fermion number arguments 
one subtracts a neutralizing background of precisely i 
charge per lattice site, whereupon {N) — > i and AA^ 
with an increasing width w. The intermediate regime 
that occurs once the zero state is covered is the crux of the 
matter. Not only does the expectation value of fermion 
number equal i, but fiuctuations are also small. After 
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the subtraction, the fermion number has the eigenvalue 
i. From the quantum optics viewpoint, this is something 
of a conjuror's trick. Correlated fluctuations in fermion 
number between adjacent sites create an impression of 
a sharp eigenvalue in a smoothly weighted sum of the 
occupation number operators for the lattice sites. 

The scattered light carries signature of the fluctuations 
in the scattered intensity. We demonstrate by plotting in 
Fig. 121 separately the contribution \E\^, as if the fermion 
numbers were precisely fixed, and the fluctuation terms 
(AA^)^. We also show the quantity {AN)f that remains 
from the fluctuation term if we only keep the contribu- 
tions with k ^ I in Eq. ((HJ, as if the fermion number 
fluctuations at adjacent sites were uncorrelated. These 
are given as functions of the width of the focus w of the 
laser beam. Here Ng = 129, Nf = 65, and we choose 
/i = 0.9 a to make a sharply localized zero state. 

Even with a very narrow focus of the laser, w — A or 
about one wavelength, the contribution from the fluctua- 
tions is two orders of magnitude below the coherent inten- 
sity, whereas the fluctuations from uncorrelated fermion 
numbers would make a contribution an order magnitude 
smaller than the coherent intensity. As our detection 
light was assumed to be far-off resonance, the photon 
number fluctuations are Poissonian. Under otherwise 
ideal conditions, the detection of about a hundred pho- 
tons could therefore reveal the difference between cor- 
related and uncorrelated fermion numbers at adjacent 
site, whereas a quantitative study of the actual corre- 
lated fermion numbers requires the detection of about 
10,000 photons. Unfortunately, a large number of scat- 
tered photons means a large number of recoil kicks on 
the fermions. Currently available optical lattices likely 
cannot absorb the assault of hundreds of photon recoils 
without developing some form of a dynamics that compli- 
cates the phenomena we are analyzing. In the coherent 
detection of the fermion number, however, the photon 
recoils could be suppressed by tuning the energy gap to 
be much larger than the photon recoil energy. 

It is instructive to note that, at the level we have dis- 
cussed (amplitude or intensity measurements), optical 
detection of the anomalously small fermion number fluc- 
tuations responsible for fractionilization has to be coher- 
ent and rely on interference of light scattered from differ- 
ent lattice sites. If a too broad angular average or other 
such cause wipes out the interferences [a^ai — s- 
we are back to adding fermion number fluctuations from 
different lattice sites as if they were independent. 

Although we, of course, do not aim at a specific exper- 
imental design, a few variations to potentially overcome 
the technical limitations wc have noted bear a mention. 
First, we have implicitly assumed that the lattice light 
and the detection light have the same wavelengths. By 
angling the beams used to make the optical lattice, or 
possibly by using microlens arrays [l7| , the optical lattice 
can be stretched. The resolution limit imposed by the 



wavelength of the detection light could be circumvented. 
Second, so far we have dealt with what in essence is spon- 
taneous Bragg scattering. Recently, induced Bragg scat- 
tering has been introduced as a method to study the 
condensates in detail P^] . How induced scattering works 
in the case of a ID lattice under inhomogeneous illumi- 
nation is not clear at the moment, but conceivably the 
Bragg pulses could be made so short that the harmful ef- 
fects of photon recoil do not have time to build up during 
the measurement. 

We have discussed optical detection of the half-fermion 
that can arise from a topological defect in an optical 
lattice holding a FD gas. Even though both the aver- 
age fermion number and its fluctuations are in principle 
amenable to optical measurements, experiments will ev- 
idently have to await further development of technology. 
In the interim, the most valuable outcome of the kind 
of an analysis we have presented would probably be the 
insights it brings into the phenomenon of a fractional 
fermion number and its prospective applications. 
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